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ABSTRACT 

Diffuse emission in the mid-infrared shows a wealth of structure, that lends itself to high-resolution 
structure analysis of the interstellar gas. A large part of the emission comes from polycyclic aromatic 
hydrocarbons, excited by nearby ultra-violet sources. Can the observed diffuse emission structure be 
interpreted as column density structure? We discuss this question with the help of a set of model 
molecular clouds bathed in the radiation field of a nearby O-star. The correlation strength between 
column density and "observed" flux density strongly depends on the absolute volume density range 
in the region. Shadowing and irradiation effects may completely alter the appearance of an object. 
Irradiation introduces additional small-scale structure and it can generate structures resembling shells 
around HH-regions in objects that do not possess any shell-like structures whatsoever. Nevertheless, 
structural information about the underlying interstellar medium can be retrieved. In the more diffuse 
regime (n(HI) < 100cm -3 ), flux density maps may be used to trace the 3D density structure of the 
cloud via density gradients. Thus, while caution definitely is in order, mid-infrared surveys such as 
GLIMPSE will provide quantitative insight into the turbulent structure of the interstellar medium. 

Subject headings: radiative transfer — turbulence — methods:numerical — ISM:dust,extinction - 
ISM:structure — infrared:ISM 



1. THE PROBLEM 

Diffuse emission in the infrared seems like a per- 
fect laboratory to study the dynamics of the interstel- 
lar medium. Recent large-scale surveys by the Spitzer 
Space Telescope, spe cifically the GLIMPSE project 
(Benjamin ct al. 2003) have provided us with unprece- 
dented high resolution data of the diffuse emission in the 
mid- infrared (MIR) . At first glance the wealth of struc- 
ture exhibited in the flux density maps seems a striking 
argument by itself for structure analysis. 

However, the conspicuous structures themselves - 
nam ely shells, bubbles, filaments and dark clouds (see 
e.g. IChurchwell et all 12004 l2006t iHeitsch et alJl2006ri 
I Jackson et all 120061 IMercer et alJ 12006(1 raise the ques- 
tion of how much of the observed structure actually corre- 
sponds to physical structure. Flux density maps contain 
information about volume density, column density and 
excitation, but to extract one of them is only possible un- 
der assumptions. For the structure analysis, ideally, we 
are interested in volume density, which is accessible only 
indirectly, leaving us with column density as a second 
best at most. In fact, over a broad range of wavelengths 
from ultraviolet to infrared, the opportunities seem to 
be rather rare where we can interpret observed intensity 
maps of diffuse emission as information about the under- 
lying column density structure. More often than not the 
medium is optically thick for the emitted radiation, or 

1 Department of Astronomy, University of Michigan, 500 Church 
St, Ann Arbor, MI 48109-1042; email: fheitsch@umich.edu 

2 University Observatory Munich, Scheinerstr. 1, 81679 Munich, 
Germany 

3 Space Science Institute, 4750 Walnut Street, Suite 205, Boul- 
der, CO 80301 

4 Department of Astronomy, University of Virginia, P.O. Box 
3818, Charlottesville, VA 22903 

5 Department of Astronomy, University of Wisconsin-Madison, 
475 N. Charter Street, Madison, WI 53706 



denser components of the ISM act as absorbers. 

Optical depth effects become weaker with increasing 
wavelength, which is why the mid-IR ( ~ 5/zm) takes a 
somewhat special position (for a study of correlation be- 
tween emission and colu mn density in the f ar-IR see e.g. 
iBethell et all 120041 and iSchnee et al.ll2006|) . At longer 
wavelengths from the far-IR through millimeter, current 
observatories have rather poor spatial resolution, pre- 
cluding study of small-scale interstellar structure. With 
a dust/PAH extinction cross section of C ex t ~ 1 0~ 23 cm 2 
per H-atom l(Drainell20"03t iLi fc Drainell2001albfl . the op- 
tical depth for MIR emission in molecular clouds should 
range below or around 1 , which would encourage a direct 
interpretation of flux density as column density. For the 
near-infrared (NI R) this possibi li ty ha s been discussed 
and supported by iPadoan et alJ i2006tl to interpret so- 
called "cloudshine" observations bv iFoster &: Goodman! 
( 2006), although there, the column densities would have 
to be substantially smaller than in the MIR. 

Although applicable to a wider range of surveys, this 
paper focuses on the MIR diffuse emission as seen by 
the IRAC camera of the Spitzer Space Telescope. To a 
large extent, the emission in the [5.8] and [8.0] bands 
comes from polycyclic aromatic hydrocarbons (PAHs) 
(to a lesser exten t they also contribute to the [3.6] 
band. iDrainel 12003). mostly excited in the environment 
of nearby UV sources. 

Two (not unrelated) issues are raised in the MIR: (1) 
The medium is generally optically thick for the soft UV- 
radiation longward of 912A that excites the PAH emis- 
sion. Thus, only the outer layers of any structures - 
unless really diffuse - are excited and thereby traced out 
in the MIR, giving the object a filamentary appearance. 
(2) Observationally, there seems to be a strong morpho- 
logical bias to shells in the PA H-emission. These coul d 
be dynamical shells (see e.g. IChurchwell et, all [2006") . 
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like e.g. windblown bubbles or HII regions. However, 
the destruction of PAHs around the UV source could 
also lead to a shell-like structure. And finally, these 
shells could be irradiation effects as noted under (1). 
And, of course, a combination of mechanisms is possible 
also. Moreover, the usual projection problem introduces 
a bias to interpreting objects as being two-dimensional, 
while they could be very "3D" : More diffuse material gets 
ionized or blown away first, leaving the ubiquitous ele- 
phant trunk remnant structures (for an impressive exam- 
ple of this problem see the GLIMPSE study of RCW49, 
lOhurchwell et~aT1l2(i?il . 

We investigate the appearance of a model molecular 
cloud in the radiation field of a nearby UV-source, in or- 
der to quantify the correspondence of various measures 
and tracers between the original column density maps 
and the derived flux density maps. Rescaling the density 
range in the models allows us to mimic different physi- 
cal environments, from diffuse clouds to dense molecular 
clouds and cores. The appearance of the original cloud 
can be completely altered by irradiation. Only for den- 
sities of up to n(HI) 100cm -3 can the flux density 
actually be interpreted as column density. Above that, 
MIR self-extinction and strong shadowing effects in the 
UV will let the maps diverge. However, even for the 
highest density range (up to n(HI) = 10 5 cm~ 3 ), where 
the flux density maps bear no resemblance to the column 
density, the structural properties of the original column 
density distribution still can be retrieved from the flux 
density. 

Our results demonstrate that the diffuse emission data 
as made available by GLIMPSE can serve as a power- 
ful means to analyze the (dynamical) structure of the 
interstellar medium. This study aims at pointing out 
possible pitfalls and at giving a rule-of-thumb estimate 
where flux density structure could be trusted to repre- 
sent column density. In the next section (<j5J), we will 
give some observational motivation in the form of diffuse 
emission maps from the GLIMPSE data. We are defer- 
ring the full structure analysis of the GLIMPSE data to 
a future paper. The models and the details of the radia- 
tive transfer treatment are described in The results 
(<JU arc summarized in SjSl 

2. OBSERVATIONAL MOTIVATION 

The following maps are extracted from images of the 
GLIMPSE project. The data were proce ssed through the 
GLIMPSE pipeline reduc tion system ijBeniamin et al.l 
l2003HWhitnev et alJ2004() . Point sources were extracted 
from each fram e using a modified version of DAOPHOT 
l)StetsonllT98 7). and the individual residual (i.e without 
point sources) frames produced by DAOPHOT were mo- 
saicked. Only sources fit well by a stellar point-spread- 
function (PSF) were extracted. Thus saturated and near- 
saturated stars remain in the images. In addition, there 
are some artifacts resulting from point-source subtrac- 
tions due to the undersampling of the stellar PSF at 
IRAC wavelengths. Figures 1-5 show [8.0] residual mo- 
saic images displayed in Galactic coordinates. We will 
split the maps into three categories: volume illumination 
(" M2 2.1|l . absorption and shadowing fi j2 2.2|) . and high- 
density environments ( W2 2.31) . 

2.1. Volume Illumination: optically thin in the UV 



Fig. 1. — Residual image of diffuse emission at [8.0] micron, 
centered on (I, b) = (51.8, +0.6). 
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Fig. 2. — Residual image of diffuse emission at [8.0] micron, 
centered on (l,b) = (343.5, -0.3). 



Figures Hand [3 show two examples of "volume-filling" 
diffuse emission, i.e. the optical-depth effect described 
in <^ does not seem to be the dominant structure gen- 
eration mechanism. Rather, it seems, the diffuse emis- 
sion traces the volume density. While Figure 2] shows 
a clear shell-like structure, Figure [21 exhibits dark fila- 
ments, probably extinction in the MIR against a bright 
diffuse background. The volume illumination could ei- 
ther stem from material that is optically thin for the 
exciting UV provided by a nearby star or by the inter- 
stellar radiation field (Fig. |2J) , or arise in an environment 
which - although at higher densities - contains several 
UV-sources, so that we still see large parts of the volume 
in PAH-emission (Fig. 

2.2. Dark Clouds: Absorption and Shadowing 

Figures |21 and 0] give two examples of structures that 
are typical of the GLIMPSE data: the combination of 
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bright and dark features in one and the same object. 

In Figured! the (Galactic) North side of the filament is 
seen in emission, while the South side shows dark struc- 
tures. Figure 0] displays a structure composed of bright 
vertical filaments, which may be an example of volume 
illumination (see previous section), while the object in 
question is located at smaller longitudes I than the fila- 
ment, at (/, b) = (309.04, -0.42): the left side is emitting, 
while the right side is seen as dark structure. In con- 
trast, the top of the same Figure shows a "pure" dark 
cloud, without any irradiation features. The reader will 
find more examples in the other frames. This raises the 
question of whether we see two physically distinct classes 
of objects here, or whether the "bright-dark" clouds are 
just an illuminated version of the classic dark clouds. 
The "bright-dark" clo uds can often be asso ciated with 
irradiating sources (see lHeitsch et aLll2006b|) . so that we 
will identify the two morphological classes as one and the 
same type of object, namely dark clouds. 

2.3. High-density Environments 

The region around (I, b) = (30.7, 0.0) (Fig. EJ shows 
one of the many bu bbles in the GLIMPSE data 
l)Churchwell et alJ2 006). Bright dust/PAH emission and 
dark clouds form an intricately structured network. The 
Northern part of the bubble exhibits several shell-like 
structures, whose brightness tapers off to the diffuse 
background more or less immediately. The region is com- 
pact and contains a lot of dark extinction regions in the 
bright emission regions, indicating that this is a high- 
density environment. Note, however, that darker regions 
in emission could also just mean that the PAHs have 
been destroyed around the UV source. Volume illumi- 
nation seems to play less of a role here, instead, the rim 
effects mentioned in the introduction seem to have come 
into full play. 

3. IRRADIATION OF A MODEL CLOUD 

3.1. The Simulations 

We employed two model series: Series A is meant to 
resemble star formation regions in the Galactic disk. It is 
derived from cloud d ispersal models in a turbulent flow 
ijHeitsch et al.ll2006a^l . and corresponds to a box length 
of 44 pc, a background density of 0.5 cm -3 , and a cloud 
density of 150 cm -3 . Model A is in approximate pressure 
equilibrium, with a temperature of 9500 K in the ambient 
medium. The model is not self-gravitating. 

Series B stems from sel f-gravitating MHD-turb ulence 
models in a periodic box ijHeitsch et alJl20~01albD . It is 
intended to be typical for conditions within a molecular 
cloud, with a box size of rs 3 pc and a mean density of 
1000 cm~ 3 , an (isothermal) temperature of 10 K and a 
turbulent rms velocity of 1 km s _1 . As we will see below, 
self-gravitation has already lead to core formation. 

Both model series A and B have a resolution of 256 3 
grid cells. The density contrast can be adapted to test 
irradiation effects in various environments. 

3.2. UV Source and Ray Tracing 

The UV-source is assumed to be a 05V star that can 
be placed anywhere in the simulation domain. Note that 
we are just investigating the radiation paths of the stellar 
UV-photons. The models do not account for ionization or 



stellar winds, and thus cannot provide us with estimates 
how they would affect the emission structure. We further 
assume that only UV-photons with A > 912A are avail- 
able for PAH-excitation, and that photons with shorter 
wavelength have already been absorbed within the HII 
region around the star. We follow the photons on rays 
through the simulation domain. The discretization in 
longitude and latit ude follows the HealPix prescription 
ijGorski et alJ l2005') in order to guarantee equal areas for 
each ray. Rays are traced within a tree-structure such 
that optical depths are known locally at each node of 
the tree. The number of rays is determined by the num- 
ber of cells on a sphere with radius r (in grid cells) as 
N ray = iirr 2 . Thus, each cell on a given sphere gets at 
least one ray. Densities are interpolated bi-linearly. 
For each location (r, 9, cf>), the energy loss is determined 

by 

dL = L(r)-L{r + dr) 

= 15 L ° ( e ~ TM " e ~ T[r+dr) ) > (!) 

with the central (stellar) luminosity Lq and the local op- 
tical depth (measured from the UV source) 

r(r) = f n(HI) C abs dr, (2) 
Jo 

where n(HI) is the atomic hydrogen particle density, and 
Cabs = 1-0 x 10 _23 cm 2 is the PAH absorption cross sec- 
tion per H-atom. The factor 4 — ' /12 accounts for the ray 
splitting: The base level consists of 12 rays, each of which 
spawns 4 children at the next refinement level /. The ra- 
dial step size is dr, thus the volume at each location 
is known, and dL can be multiplied by the correspond- 
ing volume expression, resulting in the actual UV-energy 
loss. We integrated over the whole UV-spectrum of the 
incident radiation f ield between 912 and 3000A (see e.g. 
iLi fc D raine 2 001H Fig. 1) with the constant absorption 
cross section C a bs given above. Furthermore, we assume 
that all absorbed UV-energy is converted into PAH emis- 
sio n, which we dis tribute according to the PAH-s p ectrum 
by iDraind (|200^l . As shown by ILi fc Drain j l)2001bfl 
the spectrum is independent of the radiation field over 
a large range in intensity. The emitted spectrum, sam- 
pled at 30 wavel engths, is convolve d with the IRAC re- 
sponse function ijReach et alJl2006|) for band 4 at [8.0] 
micron, which contains the bulk of the PAH emission 
l)Drainel 12003). The resulting "gray" emission is inte- 
grated along the line of sight (i.e. along a chosen grid 
axis) following the equation of radiative transfer, includ- 
ing emission and extinction, with an extinction cross sec- 
tion Cext — 1.2 x 10~ 21 cm 2 per H-atom. 

The chosen absorption and extinction cross sections 
C a bs(0.2fim) and C ex t(S/im) are a factor of sa 1 .5 smaller 
rela tive to t he va lues given by ILi fc Drainei l)2001alb| ) 
and IDraind ((2003). They give values of C a bs ~ 1-5 x 
10" 21 cm 2 and C ext 1.4 X 10~ 23 cm 2 . The slightly 
smaller cross sections can be accommodated by rescaling 
the density. Note, however, that in both cases the cross 
sections for UV and MIR differ by a factor of 100. Thus, 
material which is optically thick at UV can be optically 
thin at the re-emitted MIR. In fact, it is to a large extent 
this difference that leads to the wealth of structure in the 
observed MIR diffuse emission. 
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Fig. 4. — Residual image of diffuse emission at [8.0] micron, centered on (l,b) = (309.1,-0.4). At (I, b) = (309.00,-0.42), there is a 
"bright-dark" cloud. A "pure" dark cloud can be found at the top of the frame. 
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Fig. 5. — Diffuse emission at [8.0] micron, centered on (£, b) = (30.7,-0.0). The UV sources (sites of star formation) in the Southern 
part of the structure irradiate the surrounding remains of the molecular gas (possibly the parent molecular cloud). 



4. RESULTS 

The irradiation effects will alter the appearance of the 
cloud, depending on gas density and geometry (i j4 4.11) . 
The (point) source needs to be subtracted ( M4 4.2fl to per- 
mit a meaningful analysis of the spatial structure f ij4 4.3fl 
and to facilitate a test of the traditional structure anal- 
ysis tools, namely power spectra ( M4 4.3 4.3.1fl and struc- 
ture functions ( H44.34.3.2f) . The overall correlation be- 
tween column density and flux density deteriorates with 
increasing volume density ( H44.4|l . Under certain condi- 
tions, the flux density maps can even be used to trace 
the three-dimensional structure of the cloud ( M 4.511 . 

4.1. Morphologies 

We begin by discussing the morphology of the irradi- 
ated objects depending on the source's position and the 
density contrast. The appearance of the irradiated cloud 
depends strongly on the absolute gas density and on the 
location of the source. 

4.1.1. Irradiated Cloud 

Figures - El give a first impression of the appear- 
ance of an irradiated cloud. The top row of each figure 
shows the column density along each of the three Carte- 
sian axes. Each column stands for one of the Cartesian 
axes along which the pictures have been taken. The first 



digit of the labels in the panels gives the model series 
name (A or B). The second digit stands for the posi- 
tion of the star along the x-axis in units of L/4, where 
L = 44 pc is the box length. A value of 1 in the first 
column of images means that 1/4 of the box length is 
between the observer and the star, e.g. panel AlxO rep- 
resents a situation where the source is located between 
observer and most of the cloud, and thus the cloud "face" 
in the plane of sky as seen by the observer is irradiated. 
Thus, in the left column, the star is moved along the line 
of sight, and in the center and right columns, it is moved 
from left to right. The cloud itself is contained in ap- 
proximately the central half of the box. The third digit 
(x,y,z) gives the Cartesian coordinate axis along which 
the model has been projected. The last digit denotes the 
logarithm of the density contrast enhancement above the 
contrast already existing in the simulation, i.e. a value of 
2 denotes a density contrast enhancement of 10 2 . Thus, 
columns are a sequence in position, rows are a sequence 
in line-of-sight orientation, and each of the four Figures^ 
-[^contains frames of a fixed density contrast. The den- 
sities range between 1cm -3 at the edge of the box (in the 
"inter-cloud medium" ) and 10 2 cm~ 3 within the cloud for 
a contrast enhancement of 1, and between 1 and 10 3 cm -3 
for a contrast enhancement of 3. Question marks ("?") 
in the model names act as group specifications, e.g. A??0 
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refers to all panels in Figure EJ whereas A?xO refers to 
the panels at varying star position, projected along the 
x-axis (i.e. left column of Figure EJ) . 

For a density contrast of 1 (Fig. 0) , the cloud's column 
density structure is well traced by the PAH-cmission. 
The object is still optically thin enough for the UV pho- 
tons to excite its whole volume. Panel A3x0 and the 
row A?yO show slight self-extinction of the PAH-emission 
(slightly darkened regions within the cloud). The model 
sequence A??0 could be compared to Figures ^ and [21 

The extinction structures get more pronounced at the 
next higher density contrast of 10 or a density range 
1 < rc(HI) < 10 3 cm" 3 (Fig.0}. Panels A??l tell us more 
about these regions: The dark structures are a combina- 
tion of UV-absorption and MIR extinction. Since the 
star is moving away from the observer from panel Alxl 
to panel A3xl, the dark fuzzy region in A3xl must be lo- 
cated between the star and the observer. If it were solely 
due to UV-absorption, it should show up as a "bright" 
(light gray) region in the PAH-emission. The fact that it 
is dark in the PAH-emission maps thus shows that it is 
not only a UV-shadowing effect, but that it is also dense 
enough to extinct the background MIR radiation. A sim- 
ilar effect can be seen in panel A3zl, at the lower left of 
the cloud. The region is UV-shadowed by the bulk of the 
cloud, however, in contrast to A3xl, the optical depth in 
the MIR is smaller, so that while the background MIR 
is extincted, its brightness is close to the background 
brightness. As one can see in the corresponding column 
density map above A3zl, the region is still within the 
cloud, i.e. there remains gas available for extinction of 
the MIR. A wisp of denser gas in front of the shadowed 
region is irradiated by the source, showing up in lighter 
gray. The effect of the background intensity will be dis- 
cussed for the model series B. 

Increasing the density contrast further to 100, corre- 
sponding to a density range of 1 < n(HI) < 10 4 cm -3 , 
(Fig. [SJ panels A??2) leads to stronger absorption, and 
in some cases to complete shadowing. Simultaneously, 
the diffuse gas - still mostly optically thin for the UV 
photons - gets brighter. Overall, the maps seem to be 
more "structured" than their lower density counterparts, 
which is partly an effect of the combination of emission 
and extinction, and partly because at higher densities, 
smaller changes in density are traced out. 

That shadowing effects play a large role can be seen in 
Figures and |H1 (models A??l and A??2). Moreover, the 
spatial structure of the cloud continues from the emission 
regions to the extinction regions (panel Alz2), an effect 
often seen in the di ffuse emission maps ( see the examples 
in EOas well as lHeitsch et al.ll2006bli . 

Figure El panels A??3 stands for the extreme case: It 
corresponds to a molecular cloud with densities ranging 
up to 10 5 cm" 3 . Except for the close vicinity around 
the star, the object stays dark: all the UV is absorbed 
directly (see row A2[x,y,z]3). In panel A3?3, we actually 
catch a glimpse of the far side of the cloud: the irradi- 
ating source is on the far side of the cloud, and while 
the cloud is optically thick to the exciting UV, it is still 
(marginally) optically thin for the PAH-emission. The 
"dark fuzzy" region discussed in A3xl is still visible and 
has acquired a "halo" of excited PAHs. 

The irradiation by the central star introduces shell-like 
patterns in an object that does not have any shell struc- 



ture whatsoever (see column density panels on top of 
Figure): Especially, panel Alz3 gives the impression of a 
star sitting in a 3D-shell. In other words, it is relatively 
simple to confuse such structures with the signposts of 
interaction between a star and its surrounding medium 
(although - obviously - it does not rule out such an in- 
teraction). This situation reminds one of the shell- like 
structures North of (l,b) = (30.7,0.0), Figure^ 

FigureE3p re sents a closer look at a sub-region of panel 
Alz3 in order to learn more about the mechanisms that 
lead to the shell- like structures in the flux density maps. 
The two center panels show the column density and the 
flux density, where the central source is now located in 
the upper left corner of the frames. For an identification 
of the shell, it helps to compare this Figure to panel Alz3 
of FigureEl The top and bottom panels show column and 
flux density profiles for two rays indicated on the gray- 
scale maps. Vertical lines denote local maxima in the 
flux density, marked as well on the gray-scale maps by 
short vertical lines. Generally, we would expect emission 
peaks at locations where the radial volume density gra- 
dient is positive, and its curvature negative. However, of 
the three locations, one does not follow this expectation: 
we observe an emission peak for a negative column den- 
sity gradient. Since enhanced emission can only occur at 
positive radial density gradients, this might seem surpris- 
ing, but it only tells us that the flux density maps in fact 
trace out the 3D structure of the cloud: while the over- 
all column density drops with increasing distance from 
the source, a local density enhancement (most likely in 
the foreground because otherwise it would be masked) 
leads to the peak in the flux density. Thus, the shell- 
like structure observed here has nothing to do with a 
physical shell, but is a result of a lucky combination of 
smaller density enhancements predominantly located in 
the foreground. Observationally, such objects could be 
distinguished from physical shells by comparing them at 
other wavelengths, e.g. in the dust continuum at 24 mi- 
cron. We have to mention a caveat here, though. Our 
models do not account for PAH destruction in the vicin- 
ity of the star. This would introduce a further bias to- 
wards shell-like structures, independent of any dynamical 
effects such as stellar winds or HII regions. 

While the overall resemblance between observed flux 
density maps and the modeled ones (e.g. A??2) is sug- 
gestive, there is nevertheless one striking difference: The 
modeled maps exhibit dark radial streaks caused by the 
total absorption of UV and thus the lack of PAH emis- 
sion along those rays. Such structures are generally not 
observed, for the following reasons: (1) In the observed 
maps, the lines of sight are substantially longer than in 
the models, and thus diffuse material in the foreground 
and/or background will lead to a much higher uniform 
flux, which easily would cover any faint low-brightness 
structures. (2) The restriction to one source and neglect- 
ing the interstellar radiation field introduces a directional 
bias in our models which then leads to the shadows seen 
in the models. Only if the absorbing material is directly 
between the observer and the source can shadowing ef- 
fects be observed, namely as dark clouds. 

Clearly, the overall correspondence between column 
density maps and the PAH-emission maps deteriorates 
for increasing densities. 
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Fig. 6. — Irradiated cloud (model A??0) seen along the three coordinate axes x, y and z. The top row gives the column density in cm~ 2 . 
The star is shifted through the volume along the x-axis (see text). The flux units are arbitrary. The density contrast is 1. Shown are maps 
corresponding to the IRAC band 4. Key to the label names: First digit: model series (here: A). Second digit: stellar position in units of 
L/4 along the x-axis: A value of 1 stands for the star being located between observer and cloud (bottom left), 2 for at the box center, and 
3 for the star behind the cloud, as seen along the x-axis. Third digit: coordinate axis parallel to line-of-sight. Fourth digit: Logarithm of 
density contrast. The actual gas density ranges from 1 to 10 2 cm - . 
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Fig. 7. — Irradiated cloud (model A??l) seen along the three coordinate axes x, y and z. The top row gives the column density in cm~ 2 . 
The star is shifted through the volume along the rr-axis (see text). The flux units are arbitrary. The density contrast is 10. Shown are 
maps corresponding to the IRAC band 4. Key to the label names: First digit: model series (here: A). Second digit: stellar position in units 
of L/4 along the x-axis: A value of 1 stands for the star being located between observer and cloud (bottom left), 2 for at the box center, 
and 3 for the star behind the cloud, as seen along the x-axis. Third digit: coordinate axis parallel to line-of-sight. Fourth digit: Logarithm 
of density contrast. The actual gas density ranges from 1 to 10 3 cm - . 
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Fig. 8. — Irradiated cloud (model A??2) seen along the three coordinate axes x, y and z. The top row gives the column density in cm~ 2 . 
The star is shifted through the volume along the x-axis (see text). The flux units are arbitrary. The density contrast is 100. Shown are 
maps corresponding to the IRAC band 4. Key to the label names: First digit: model series (here: A). Second digit: stellar position in 
units of L/4 along the x-axis: A value of 1 stands for the star being located between observer and cloud (bottom left), 2 for at the box 
center, and 3 for the star sitting behind the cloud, as seen along the x-axis. Third digit: coordinate axis parallel to linc-of-sight. Fourth 
digit: Logarithm of density contrast. The actual gas density ranges from 1 to 10 4 cm~ 3 . 
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Fig. 9. — Irradiated cloud (model A??3) seen along the three coordinate axes x, y and z. The top row gives the column density in cm~ 2 . 
The star is shifted through the volume along the x-axis (see text). The flux units are arbitrary. The density contrast is 1000. Shown are 
maps corresponding to the IRAC band 4. Key to the label names: First digit: model series (here: A). Second digit: stellar position in units 
of L/4 along the x-axis: A value of 1 stands for the star being located between observer and cloud (bottom left), 2 for at the box center, 
and 3 for the star behind the cloud, as seen along the x-axis. Third digit: coordinate axis parallel to line-of-sight. Fourth digit: Logarithm 
of density contrast. The actual gas density ranges from 1 to 10 5 cm - . 
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Fig. 10. — Model Alz3 in a zoomed-in view. The central source 
is now located in the upper left corner. Over-plotted are two rays 
through the column density and flux density map respectively, 
whose profiles are shown in the top and bottom panel. Vertical 
lines denote points of local maxima in the flux density along the 
ray. In order to help identification in the maps, these points are 
marked with small vertical lines on the ray. Not all flux density 
maxima have corresponding column density maxima. To identify 
the shell-like structures more easily, compare to panel Alz3 in Fig- 
ure ini 



4.1.2. Embedded Star 

Figure ITT1 is a variation on the previous theme, namely 
flux density maps of a self-gravitating turbulent medium, 
irradiated by a star. As before, a sequence in increasing 
density contrast is shown. One difference between model 
A and B is obvious: Model B does not exhibit a well- 
defined irradiated cloud structure (or at least the rims of 
that). One could see model B as a "close-up" of model 
A, taken of a small region around the star inside the 
cloud of model A. In model B, the UV-radiation is mainly 
stopped by the numerous filaments arising from shock 
compressions. In addition to those filamentary struc- 
tures, the emission has a "fog-like" appearance: in con- 
trast to model A, the density within the cloud is already 
high enough to absorb sufficient UV to trace even low- 
density regions, and the scales of the frame are smaller 
by a factor of nearly 15, so that we are seeing the close 
environment around the O-star. The small black objects 
visible at the higher density contrasts are self-gravitating 
cores that have already collapsed. While it is possible to 
identify irradiated filaments with structure in the col- 
umn density map (top), the overall agreement between 
flux density and column density (even for the lowest con- 
trasts) is rather poor. The densities range between 10 
and 10 4 cm~ 3 in the original density distribution (density 
contrast enhancement of 0, model B??0). An enhance- 
ment of 1 then corresponds to a maximum density of 
3 x 10 4 cm- 3 , 2 to 10 5 cm~ 3 and finally 3 to 3 x 10 5 cm- 3 . 

The dark regions in the lower row of Figure 1111 are 
reminiscent of the dark clouds in Figure however, for 
the wrong reason: the observed dark regions arise from 
back-illumination by the diffuse 8 micron emission along 
the whole line of sight, while the dark ray-like regions in 
model B are a result of the complete UV-absorption of 
the single source available in the model. 

Model sequence B does not seem representative of ob- 
served diffuse MIR emission. The column density maps 
prominently exhibit filamentary structure over the whole 
domain. These filaments are shock fronts caused by the 
supersonic turbulence in the box. However, they are not 
reproduced in the flux density map. Obviously, because 
of their highly transient nature, the filaments contain 
only little mass, so that they do not represent a major 
obstacle to the star's UV-radiation. Stars form in these 
models predominantly in the (abundant) intersections of 
filaments. The high shock compression leads to a close to 
instantaneous destabilization and fragmentation of these 
intersections, resulting in the formation of collapsed cores 
afte r less than one dynamical time in the simu lations (see 
e.g. iKlessen et alJl2000l iHeitsch et al.ll20f)laj) . In other 
words, cores form too soon in these simulations for the 
filaments to gather sufficient mass to be seen in MIR 
diffuse emission. 

4.2. Intensity Distributions 

Before we can analyze the intensity distributions (this 
section) and the spatial structure f H4 4.3Jl , we need to "re- 
move" the central source, which would otherwise domi- 
nate the analysis. To that purpose, we average the inten- 
sities azimuthally, centered on the flux density peak, and 
subtract the resulting profile from the frame. An alter- 
native approach would be to just cut off the peak of the 
intensity corresponding to the central source. However, 
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Fig. 11. — Irradiated cloud (model B) seen along the three coordinate axes x, y and z. The top panels of the upper row give the column 
density in cm -2 . The star is shifted through the volume along the a;-axis (see text). The flux units are arbitrary. From top left to bottom 
right, the density contrasts are 1, 3, 10 and 30. Shown are maps corresponding to the IRAC band 4. Key to the label names: First digit: 
model series (here: B). Second digit: position in units of L/4. Third digit: coordinate axis parallel to line-of-sight. Fourth digit: Logarithm 
of density contrast. 



this leaves the imprint of the i? _2 -dilution of the radia- 
tion field, which is a direct effect of the central source and 
would confuse the structure analysis. The subtraction of 
the i? _2 -profile could be seen as a modified photometry 
algorithm for point sources irradiating a surrounding, op- 
tically thin diffuse medium. 
The residual frame (Fig. I12|) still shows artifacts due 



to the central source, but the cloud structure is more 
pronounced than in the "directly observed" maps. Com- 
paring the residual maps to the column density maps 
(top row of Fig. EHHI) demonstrates that the irradiation 
- besides introducing additional structure via absorp- 
tion (dark regions) - generates more small-scale struc- 
ture than observable in the column density map. This is 
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Fig. 12. — Residual flux density maps for the three lower density 
contrasts of model A. Artifacts unfortunately are unavoidable. 



mostly an optical depth effect in the sense that already 
small increases in the optical depth (shocks, filaments, 
smooth gradients in the density) leads to additional ab- 
sorption which then will be traced out by the re-emitted 
MIR. This effect gets amplified if the "obstacle" uses up 
all the remaining UV photons in that ray, i.e. if the op- 
tical depth changes from t < 1 to t > 1 within that 
structure. The result is a bright region facing the central 
source, and a dark, shadowed region facing away from 
the source. 

Subtracting the central source in model B (Fig. H3jl 
is not as prone to produce artifacts in the residual flux 
density maps because the model is pretty much isotropic 
by construction. In contrast, since the cloud center of 
model A does not coincide with the UV-source, artifacts 
are unavoidable and easily recognized as darker circular 
regions. 

4.3. Spatial Structure 

In the previous section we saw that a "by eye"- 
comparison of the column density maps and the flux den- 
sity maps reveals additional structure in the flux density. 
The next two sections will discuss under which condi- 
tions, and how strongly, the irradiation affects structure 
measures. We test the reliability of two traditional struc- 
ture analysis tools, namely power spectra ( H44.34.3.1|l 
and structure functions ( ^4 4.3 4. 3. 2J) . Since the flux den- 
sity maps at higher densities are dominated by global 
irradiation effects, power spectra turn out to retrieve 
structure information less reliably than structure func- 
tions. The latter can only be trusted if applied to small 
enough regions not to be influenced by global irradiation 
effects. 

4.3.1. Power spectra 

How reliable is the structure information extracted 
from diffuse emission? Figure 1141 (left column) gives an 
impression for our irradiated cloud (model A). The top 




Fig. 13. — Residual flux density maps for the three lower density 
contrasts of model B. Since the model tends to be more isotropic 
by construction, artifacts of the subtraction carry less weight. 



panel shows the power spectrum for the original column 
density map, without any source or irradiation. The 
three lower panels in turn give the spectra of the irra- 
diated cloud including the source (solid lines on top of 
panels), and the spectra of the residual flux density, i.e. 
with the central source subtracted as described in M4 4.21 
Obviously, the maps including the source are completely 
useless for structure analysis: the source dominates the 
whole spectrum. The slope ranges around —0.2 (not 
shown in the plot). A completely different picture is re- 
vealed when analyzing the residual maps: The resulting 
slopes approach the values of the column density maps. 
The slightly flatter slopes indicate that irradiation in- 
troduces additional small-scale structure. The artifacts 
introduced by the subtraction of the central source in 
model A2?0 are mirrored in the slight hump of the cor- 
responding spectrum around logfc « 0.9. Note that the 
power law exponents refer to power spectra, i.e. a Kol- 
mogorov law would be represented by —10/3. 

To conclude for model sequence A, the power spectra 
reproduce the original scale distribution within the er- 
rors, which however are substantial, due to the circular 
averaging. 

The situation changes drastically for model series B 
(Fig. El right column), because the high-density cores 
lead to a strong power increase at small scales (basically, 
they act as noise) in the column density spectrum, which 
is why this spectrum is much flatter than expected for a 
turbulent power spectrum (see also lHeitsch et al"ll2001al 
for a discussion of this effect). Note that the column den- 
sities in the cores range approximately 2 orders of mag- 
nitude above the mean column density (Fig. ITTJ) . Since 
only a tiny fraction of a dense core can emit in the MIR, 
these objects generally do not show up in the emission 
maps, and then only in extinction. Again, the central 
source dominates the "uncleaned" spectrum. Thus, in a 
sense, the MIR maps of model B are more suitable for 
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Fig. 14. — Power spectra of column density (top) and flux density (remaining panels) for models A2?0 through A2?2 (left column) and 
B2?0 through B2?2 (right column). All spectra are normalized to the mean. Diamonds stand for spectra of the residual maps, while the 
solid lines around logP(fc) = correspond to the spectra of the maps including the central source. The fitted power law exponents a XyVyZ 
are given in the panels, where the indices denote the direction of the line-of-sight (e.g. a x belongs to A?x2). The error bars denote the 
logarithmic error on the mean, while the errors in a contain the fit error as well as the error on the mean. 



structure analysis than the original column density maps 
- granted that the major part of the map is unaffected 
by shadowing or self-extinction. 

4.3.2. Structure Functions 

Structure functions are a more appropriate tool for an- 
alyzing the spatial distribution of an incomplete data set, 
since masking introduces artificial signals in the power 
spectra. Only pixels with a measurable signal would 
contribute to structure functions of observational data. 
However, due to the irradiation and absorption effects 
the flux density map generally will not cover the whole 
cloud, as demonstrated in H44.ll Thus, in order to com- 
pare the column density and flux density maps, we se- 
lected the former for column densities N > 10 22 cm -2 , 
which traces out the bulk of the cloud. Applying the re- 
sulting mask to the flux density maps and determining 
the structure functions 

SF(l) = (\F(x)-F(x + l)\ 2 ) x (3) 

with lag I and spatial coordinate x, results in Figure IT5l 
The four top rows belong to model Al??, with the 
source being offset to smaller x. The three center rows 
refer to model A2??, with the source dead-center, and 
A3?? is represented by the four bottom rows, where 
the source is offset to larger x, so that A3x? refers to a 
configuration where the cloud is located between observer 
and source. Each panel shows the structure function of 



the original column density at N > 10 22 (thin line). The 
structure function for the corresponding flux density is 
plotted in thick line style. The numbers s in each panel 
refer to the logarithmic slope of the structure function, 

logSF(l)<xl s (4) 
taken over the first decade in I. Although a column den- 
sity threshold of N > 10 22 leads to most of the high- 
density gas in one contiguous object, there are a few iso- 
lated cloud fragments. Since these cloud fragments are 
separated from the main body of the cloud, they tend to 
have the largest lags with respect to the main body, how- 
ever, it turns out that their column densities are close to 
the mean column density within the main body, so that 
the structure functions show a drop at large lags. Since 
Figure lT51 is somewhat messy to interpret, Figure H^l sum- 
marizes the degree to which the original column density 
structure information is retrieved from the flux-density 
maps 

Just looking at the filled symbols (averages over the 
three coordinate directions) one might get the impression 
that larger density contrasts lead to more accurately re- 
trieved slopes. This is somewhat surprising. However, 
the scatter between the three coordinate directions is 
substantial and suggests to take a closer look at Fig- 
ures EHH1 again. 

As will be shown in H4 4.5I the lowest density contrast 
allows a "volume-rendering" of the cloud: the medium is 
optically thin enough so that the exciting UV can light 
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Fig. 15. — Structure functions of column density N (thin lines) 
and residual flux density maps (thick lines) for the three spatial 
directions. The center block shows models A2??, i.e. with the 
source dead-center, while the upper and lower block correspond to 
Al?? and A3?? respectively, i.e. the source offset to smaller/larger 
x. Since model A2?3 would not be very revealing (see Fig.|2j, wc 
dropped it. 



up the whole cloud, nearly acting as a uniform back- 
ground field. If the source is within the cloud, though, 
it tends to over-represent the volume directly around it. 
Thus, the best fit is given by A3y0, where the source 
is located seemingly at the tip of a filament, farthest 
away from the bulk of the cloud. The effect of over- 
emphasizing close-by structures can be well observed 
even in A3y0: above the source, a filament is traced 
out which does not show up clearly in the column den- 
sity map. Comparing this to A3x0 we see that it suffers 
from "neighborhood-illumination" by the source: inter- 
nal structure of the cloud thus is just smeared out: the 
flux-density slope is flatter than the column density slope 
(Fig. I15fl . A3z0 fares slightly better, although here the 
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Fig. 16. — Normalized difference between logarithmic slope of 
structure functions (eq.|3} for column density (sjv) and flux density 
(sp), against the density contrast. For reasons of clarity, the three 
coordinate directions are plotted in open symbols centered around 
their respective density contrast, i.e. A?xl can be found left of 
An/Ano = 1, A?yl at 1, and A?zl right of 1. Filled symbols 
denote the mean over the axis directions x, y and z. Symbol types 
refer to the position of the source. 



source over-emphasizes the ring-like structure close by, 
leading to a steeper slope. The center-position maps 
(A2?0) all suffer from the source being located inside the 
cloud, while models A1?0, where the source (along the 
a;-axis is located in the front, again tell a mixed story. 
For AlzO, the source is located at the upper left end 
of the cloud, in a fairly diffuse environment and suffi- 
ciently far away from the bulk of the cloud in order not 
to affect the structure function (third cross symbol for 
logAn/Ano = in Fig. fT5|l. For AlxO and AlyO, the 
source is already too embedded. 

Moving on to the next density contrast, log An/Ano = 
1, again, A3yl shows the lowest deviation, however, be- 
cause of the shadowing effects in the cloud's left ("East- 
ern") part, the deviation is larger than for A3y0. That 
A3xl provides better fits than A2xl and Alxl is under- 
standable again because of shadowing, although it is not 
directly clear why the match should be better than that 
for A3x0. Checking the corresponding panels A3x0 and 
A3xl in Figure IT51 gives the answer: on the small scales, 
the structure functions agree, however, on large scales 
they differ. The better matches for Alzl and A2zl (com- 
pared to A3zl) might be a consequence of the shadowed 
regions, which still show substructure: this substructure 
would enter the structure function on the small scales 
and provide the same information as structure in emis- 
sion. Thus - granted that we chose a small enough field 
for analysis - to some extent structure in absorption and 
emission can be used simultaneously for analysis. 

At density contrasts log An/Ano = 2 and 3, match- 
ing between the overall shapes of structure functions for 
column and flux density deteriorates on the large scales 
(Fig. HH|. as one would expect by looking at Figures 00 
and[!|] However, there is still enough signal from irradi- 
ating small density variations for the structure functions 
to match on the small scales. Again, given that we take 
a small enough region, the structure function can be re- 
produced. 

Another limitation on the size of the region is given by 
the available dynamic range in the observed flux density 
maps. In the diffuse emission, it often turns out that the 
dynamic range is spanning 1, at most 2, decades, which 



16 



Structure Generation by Irradiation 





0.6 


■ ft 
















'. ft 






ft 

o 


A??2 
B??2 




- 




0.4 
















'(N,F)> 


0.2 


" & 












- 


S3 

V 




: o 
















0.0 
-0.2 


- o 
o 


o- 


o 

# 









ft 


ft ; 






0.0 


0.5 


1.0 


1.5 




2.0 


2.5 3.0 



log An/Ano 

Fig. 17. — Correlation coefficient (eq. |FJ for model series A and 
B. For larger density contrasts, the maps are uncorrelated. 



limits the spatial scales available for analysis. 

To summarize, structure functions seem to be a much 
more viable tool to investigate localized diffuse emission 
than power spectra. Despite the fact that the irradiation 
sometimes modifies the underlying density information 
beyond the point of recognizability, the resulting struc- 
ture functions still retrieve the salient scale information - 
given that the field investigated is small enough. Because 
of the global irradiation effects, there is only little hope to 
retrieve the large-scale information accurately. The most 
promising way would be to identify small, "uncontami- 
nated" regions and reconstruct the structure information 
from those. 

The deviations (Fig. I16fl are within the errors on the 
mean of the structure function (see the analysis of power 
spectra). Again, these stem from the circular averaging. 
The large deviations in fact point to a major problem 
when applying two-point correlators to data exhibiting 
anisotropic structures. For example, filaments will lead 
to different structure information along and perpendic- 
ular to their long axis, in turn leading to a mixing of 
the spectral information when averaging the directions. 
Since the observed structure in the interstellar medium 
seems to be anisotropic to a large extent, the application 
of unmodified two-point correlators might raise ques- 
tions. We will discuss this problem and possible remedies 
in a future paper. 

4.4. Correlation Measures 

How reliably are structures in column density repro- 
duced in the flux density maps? The correlation coeffi- 
cient 



(C(N,F)) = 



(5) 



serves as a first crude estimator for the similarity of col- 
umn density maps N{j and flux density maps Fij. The 
means over the map are denoted by (N) and (F) respec- 
tively. For all density contrasts larger than 1, N and F 
are basically uncorrelated fFig. lTTjl . 

The correlation coefficient provides an overall estimate 
for how "similar" the maps are. Unfortunately, equa- 
tion Q is only of limited use for a point-to-point sim- 
ilarity measure. Although one could plot a map of the 
summands, the normalization is unclear. One possible 




Fig. 18. — Correlation maps (eq. |HJ for model series A. 



similarity measure is given by 

N — (N) F - 



C(N,F) = 



(F) 



(6) 



(N) (F) ' 

At locations where both maps show either strong or weak 
signals, C > 0, while positions with anti-correlated sig- 
nals will exhibit C < 0. 

Figure ED shows C(N, F) for model A. It can be di- 
rectly compared to Figures H3IE1 and El For the lower 
density contrasts, the column density map correlates well 
with the emission (light gray to white shades), while for 
higher density contrasts, the absorption regions lead to a 
strong anti-correlation. The overall positive correlation 
deteriorates for larger density contrasts. 

4.5. PAH-Emission as Edge Sensor 

If the UV from the central source excites PAH-emission 
in the rims of the surrounding molecular cloud, the emit- 
ted flux density could be used as an edge sensor or gra- 
dient indicator. To test this, for each position in the 
simulation domain and the column density map we cal- 
culate the 3D radial density gradients 



P- P 
AR 



(7) 



and the two-dimensional column density gradients re- 
spectively 



d r N(rij,<l>ij) 



N - N 
Ar 



(8) 



1,3 



The hatted quantities are bi-linear interpolations of the 
(column) density at the positions 

p= p(Xi +AR sin 9 cos <f), 

Yi +AR sin 9 sin cf), 

Z l +ARcos9) (9) 

N = N(xi +Ar COS0, 

j/ 4 +Arsin0), (10) 
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with R = \/X 2 + Y 2 + Z 2 and r = ^Jx 2 + y 2 . The ra- 
dial increments AR and Ar measure one grid cell, and 
6 as well as 4> m equations 10 and (|ll)fl refer to the co- 
ordinates (i,j,k), (i,j) respectively. The gradients are 
computed with respect to the central source (R = and 
r = at position of source) and selected for positive val- 
ues (Figs 120)1 ■ Positive gradients dnp > indicate 
locations of enhanced UV-absorption and thus enhanced 
PAH-emission, while d r N would denote the same if the 
physical cloud structure were truly 2D. We compute dap 
and d r N for each cell in the simulation domain and col- 
umn density map respectively, i.e. A slight subtlety en- 
ters regarding the 3D-case: We still need to project the 
gradients on a plane. In order to keep the 3D-information 
at least partially, we color-coded the distance along the 
line of sight, where blue colors denote short distances 
(towards the observer) and red colors denote long dis- 
tances (at the far end of the volume). The intensity 
gives the strength of the gradient. The distance infor- 
mation is weighted with the gradient strength, i.e. two 
gradients of same strength located at one quarter and at 
three quarters of the box length will show up in purple. 

At first sight, the correspondence between gradient 
maps and flux density is at most poor. A close look 
however reveals some interesting details. We focus on 
Figure HI (models A2?l and A2?2) first. The top row 
shows the original column density, followed by the 2D 
gradient d r N and the projected 3D gradient dnp. The 
two bottom rows give the flux density maps for mod- 
els A2?l and A2?2, i.e. for increasing density contrast. 
The color scale is a measure of the location of the gra- 
dient, with blue closest to the observer, and red at the 
far end of the box. Clearly, the gradients making up the 
map come from all possible locations along the line of 
sight. Taking A2yl, the L-shaped structure visible in 
extinction in the flux-density maps is traced out by a 
hazy bluish gradient-indicator: The corresponding den- 
sity jump is close to the observer, which makes sense, 
since dense material between the source and the observer 
leads to the extinction. Further right of the L-shape, the 
gradient map of A2y? shows a red (albeit dark) region: 
weak gradients at the far side of the source. Correspond- 
ingly, we see a bright irradiated region in the flux den- 
sity map of A2yl: the (far away) structures are irradi- 
ated by the central source. This is not seen anymore in 
A2y2: the density has increased and leads to extinction 
of the far-away emission. Similarly, the dark region to 
the lower left of the center in A2xl corresponds to gra- 
dients close to the observer. In this case however, there 
is no clear blue signal in the 3D gradient map, which in- 
dicates that the dark region is caused by a combination 
of UV-shadowing and MIR extinction of diffuse back- 
ground. Moving northwards of that direction, the source 
of emission is receding from the observer: again, this in- 
formation is lost at higher density contrast (A2x2), and 
thus the cloud structure is only incompletely "sampled" . 
Interestingly, for A2z2, the situation is nearly converse: 
The dark gradient-regions link up with strong emission 
signals (to upper left of the center), while the extincted 
regions to the South are farther away from the observer. 
Thus, those dark regions might again be caused by a 
combination of missing irradiation and extinction of dif- 
fuse background. The "handle" to the right is located at 
the far side of the central source. 
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Fig. 19. — Gradient maps (see text) for model A, source centered 
(A2??). Column density (top row) and flux density (two bottom 
rows) are shown for comparison. The color bar for the 3D gradient 
gives the distance along the line of sight in box units, with 0.0 
(blue) closest to the observer, and 1 (red) located at the far end of 
the box. 



Compared with this richness of structure, the 2D gradi- 
ents come as sort of a big disappointment: They resemble 
the cloud structure only in the widest sense. Beginning 
with A2yl again, the L-shape is recovered, but this is al- 
ready all we can see. For A2z?, the "handle" to the right 
of the central source is showing up slightly in d r N. How- 
ever, moving to the larger density contrast (A2?2), the 
situation changes slightly: The flux density maps look 
more "rim" -like, more easily identified with the 2D gra- 
dient map. Carrying this idea to the extreme, Figure I2TH 
shows the gradient maps and flux densities for Al?3 (as 
Fig. shows, model A2?3 would be not exactly informa- 
tive in this context). The 3D gradients now are hardly 
reproduced (with the noticeable exception of the South- 
ern absorption rim in Alx3). The 2D gradients however 
tend to trace now the absorption rims (see e.g. Aly3 
and Alz3). Thus, for lower density contrasts, or a more 
diffuse environment, the emission structures are intrin- 
sically 3D, while for higher density regions, they tend to 
get more and more 2D (although 3D effects are still not 
to be neglected) . 
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Fig. 20. — Gradient maps (see text) for model A, source at 
position 1 (Al??). 



5. SUMMARY 

The abundance of structure in MIR diffuse emission as 
observed in e.g. the GLIMPSE data seems to offer a per- 
fect laboratory to study the dynamics of the dense ISM. 
However, a large part of the observed structure could be 
irradiation effects due to PAH emission. PAHs are ex- 
cited by UV photons from nearby stellar sources or the 
interstellar radiation field, and re-emit in the MIR. Be- 
cause the respective cross sections differ by a factor of 
approximately 100, gas, which is optically thick in the 
UV, can still be optically thin in the (emitted) MIR. 
Thus, PAH emission is often seen in filamentary struc- 
tures, probably the "rims" of denser clouds. These ir- 
radiation effects might spoil the opportunity to study 
the ISM gas structure, since this requires interpretation 
of the observed flux density in terms of volume or col- 
umn density. Motivated by a few examples taken from 
GLIMPSE, we identified possible limitations of this in- 
terpretation. We quantified the reliability of flux-density 
maps of diffuse emission in the MIR to reproduce the 
underlying (column) density information. We used two 
model sets, one corresponding to a (more or less) isolated 
"molecular" cloud, and the other imitating a region deep 
inside a molecular cloud, both irradiated by an 05 star. 
PAHs absorb the UV and re-emit the energy in the MIR, 
which then is integrated along the line-of-sight, including 
MIR extinction. In the following, we will summarize how 
reliable flux density maps reproduce column density in 
our models, and how this affects the structure analysis, 
with possible applications to MIR observations. 

5.1. Morphology: column density and flux density 

If the medium is optically thin for the irradiating UV, 
then the MIR emission maps could be used for a high- 
resolutio n study of the colum n density structure of the 
medium. iPadoan et alJ l)2006l) supported this possibility 
fo r the NIR, based on the o bservations of "cloudshinc" 
bv lFoster fc Goodmanl (2006) . Since they were interested 



in the appearance of the diffuse emission, they used an 
isotropic radiation field mimicking a UV background, in 
contrast to our models that employ a point source for 
irradiating the surrounding medium. For higher-density 
environments such as molecular cloud in the vicinity of 
strong UV sources, interpreting MIR emission as column 
density requires some caution: The transition from r < 1 
to r > 1 can lead to strong signals in the MIR flux den- 
sity, but not necessarily paired with a corresponding sig- 
nal in the column density: The maps bear no resemblance 
to the column density (Figs IBH^l and ITTJl . In the extreme 
case, PAH-emission is only excited at the rims of clouds 
(e.g. iChurchwell etaTfe OOfil. causing the impression of a 
highly filigree structure in the gas: irradiation introduces 
more small-scale structure than observable in the under- 
lying column density maps. As soon as shadowing is 
obvious, the structure seen in emission will generally not 
represent column density. Some of the observed shell- 
like structures could be just irradiation effects, and by 
themselves indicate that shadowing (or "rimming") has 
set in (Figs. 151 171 and fT0j> . Since our models do not in- 
clude PAH destruction around the star, we expect them 
to exhibit fewer shell-like structures than observed. PAH 
cavities would lead to "shells" even if the cavities were 
not associated with e.g. wind-blown bubbles. 

5.2. Structural Properties 

Power spectra are only partially useful for an anal- 
ysis of diffuse emission structure. Their well-known 
main drawback is that they tend to confuse the informa- 
tion about extent of a region and the separation of re- 
gions. Furthermore, masking is always an issue in power 
spectra, since it tends to introduce a signal by itself. 
Power spectra containing the central source are com- 
pletely dominated by that fFig. 1141 left column), while 
spectra of residual maps are slightly flatter than the un- 
derlying column density distribution. This could be a 
projection effect and/or the result of additional small- 
scale structure traced out by irradiation. The spectral 
slope is pretty much insensitive to the density contrast 
within the error bars, which are significant. For collapsed 
regions fFig. 1141 right column), the column density spec- 
trum flattens considerably because of the strong point 
source contribution. Compared to that, the flux density 
spectra steepen because the point sources are not fully 
irradiated and thus do not show up (except in extinc- 
tion) . 

Structure functions seem to be a more viable tool to in- 
vestigate localized diffuse emission. Despite the fact that 
the irradiation may modify the underlying density infor- 
mation beyond the point of recognizability, the resulting 
structure functions still retrieve the salient scale informa- 
tion - given that the field investigated is small enough not 
to be contaminated by global irradiation effects. There 
is little hope to retrieve the large-scale information ac- 
curately by applying a global structure measure such as 
power spectra. Structure in extinction can be used as a 
continuation of structure seen in emission, although this 
raises the issue of an appropriate choice of background 
for the extincted region (Figs. and [7J. The deviations 
(Fig.llfi| are within the errors on the mean of the struc- 
ture function. 

Since the overall structure in the ISM tends to be 
anisotropic, applying two-point correlators seems at least 
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questionable. Averaging in k or I space leads to substan- 
tial errors on the mean, which themselves indicate that 
the underlying structure is anisotropic to a large extent. 

The application of these results to actual observational 
data (GLIMPSE) and the discussion of anisotropy we 
defer to a future paper. 

5.3. Flux density as gradient indicator 

Since the conversion of UV to MIR will occur predom- 
inantly at regions of large positive radial density gradi- 
ents (as long as there are photons left), the flux density 
maps might offer the opportunity to gather information 
about the (3D) density structure of the cloud. To test 
this, we compared the radial volume and column density 
gradients to the flux density maps. For lower density 
contrasts, the flux density maps tend to trace out the 
3D structure of the cloud, and in fact they can be used 
as 3D gradient indicators. For higher density contrasts, 
they revert to an indicator of the column density gradi- 
ents: the structures seemingly become two-dimensional. 
Thus, more diffuse regions are intrinsically "more 3D", 
while higher-density environments tend to be 2D. 

5.4. Limitations 

Beside the irradiation effects discussed here, there are 
other limitations to an interpretation of flux density 
maps as column density. 

(1) If the volume density is low enough that the excit- 
ing UV can irradiate the whole cloud, one might question 
how long the line-of-sight actually is, and whether angu- 
lar effects leading to scale-mixing in a structure analysis 
would play a role. 

(2) On the other hand, the volume density might be 
large so that the irradiating UV will be absorbed more 
or less directly "at the rim" (if such a thing exists) of the 
cloud. Then, depending on the geometry of observer, 
irradiated medium and irradiation source, the observer 
might see predominantly ID structures or 2D structures, 
implying projection effects in the spectral information. 

5.5. Conclusions 
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Depending on the diagnostics, MIR flux density maps 
of diffuse emission from PAHs excited by a nearby UV 
source can be used to extract information about the 
density structure of the underlying (molecular) cloud, 
though this statement needs some qualification. 

(1) Flux density maps need not correspond "by eye" 
to column density maps: due to irradiation effects they 
tend to show more small-scale structure. 

(2) Irradiation by a point source can produce shell- 
like structures, mimicking physical shells, even in objects 
which do not have any shell-like properties. 

(3) As long as structure studies are restricted to ar- 
eas small enough not to be contaminated by any large- 
scale effects, flux density and column density show sim- 
ilar structural properties. However, the application of 
unmodified two-point correlation functions introduces 
substantial errors on the mean due to the underlying 
anisotropy in the ISM structure. 

(4) MIR flux density maps tend to trace out gradients 
in the three-dimensional density distribution. 

Employing MIR diffuse emission to extract structure 
information about the underlying interstellar medium re- 
quires close attention to the environment. This study at- 
tempts to provide some guidelines to chose appropriate 
locations. Bearing the limitations in mind, analyzing the 
ISM structure with the help of the GLIMPSE data will 
be a promising task. 
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